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Abstract 

We report on a measurement of the decay Kl /x"'"/x~7 from Fermilab 

experiment E799. We observe 207 candidate signal events with an estimated 

background of 10.5 ± 4.0 events and establish B{^Kl /x"'"/x~7) = (3.23 ± 

0.23[stat) ± 0.19[sys)) X 10~^. This provides the first measurement of the 

Kjj* form factor in the muonic Dalitz decay mode of the K^. 
PACS numbers: 13.20Eb, 14.40Aq. 

Accepted for publication in Physical Review Letters. 



2 



Typeset using REVTeX 



The decay fi'^fi'^f is expected to be dominated by the so-called long distance 

contributions occurring via the K"f"f* vertex [1]. Knowledge of long distance contributions 
can help to determine the origin of the AI = | enhancement in non-leptonic weak inter- 
actions. Other rare kaon decays such as depend on similar long distance 
contributions which must be measured in order to extract standard model parameters [2-4]. 
The first study of the decay Kl n'^ is an important test of the current theoretical 
understanding of long distance effects in neutral kaon decays. 

The differential decay rate for Kl l'^l~"f can be calculated using QED modified by a 
form factor f{x) for the K"f"f* vertex. The result normalized to Kl "f"f is given by 

"ax Sir X 



1-xr 1 + ^ 1 



xra\j \ xra\ ^ 



where x = mfi/mj^, and the normalization is such that /(O) = 1. Bergstrom, Masso, 
and Singer (BMS) [1] parametrized this form factor using a combination of a vector meson 
dominance model [Kl ir°,r],r]' "f*"f) and a vector- vector transition model [Kl 
K*"f [p,io, (j))*"f "f*"f)- In this model there is a single free parameter ax* and f[x) = 
fi[x) + O-K* fK[x), where fi[x) represents the vector meson dominance amplitudes and fx 
is a sum over pole terms from a model of the KK*"f vertex. The parameter ax* has been 
inferred from studies of the decay mode Kl e"'"e~7 [5,6]. The model of Ko [4] has no 
free parameters, and makes a definite prediction for the form factor. At this point chiral 
perturbation theory gives no prediction for Kl n'^ as diagrams cancel at 0[p^) and 
higher order terms have not been calculated [7]. 

In the decay Kl l'^l~"f the di-lepton invariant mass distribution is constrained to be 
between 2m/ and mx- Together with the 1/x dependence of the differential rate in Eq. 1 
this gives the electron mode a much larger region of decay phase space as compared to the 
muon mode. Thus measurements of the form factor f[x) will be most sensitive to the di- 
electron mass spectrum (and not the Kl e~^e~"f branching fraction) and the Kl n'^ 
branching fraction (but not the Kl n'^ di-muon mass spectrum) assuming uniform 



experimental acceptances. 

Fermilab experiment E799 was designed to search for rare CP violating decays in- 
cluding Kl 7r°e"'"e~ and tt^ n'^ n~ [8]. In the fixed target mode of the Tevatron, 
800 GeV protons were incident on a beryllium target. Two nearly parallel neutral beams 
were produced including if^'s which then decayed in fiight. The detector is described in 
detail elsewhere [9,10]; only those elements relevant to this analysis are mentioned here. 
Two hodoscope banks of vertical and horizontal counters 175 m downstream of the target 
were used to trigger on kaon decays with two or more charged tracks. Just downstream of 
these was a circular calorimeter of 804 lead-glass blocks with cross sections 5.8cm X 5.8cm 
and 18.7 radiation lengths long. The resolution for typical photon energies (12.5 GeV) in 
the calorimeter was 5%. Four drift chambers each with single hit resolutions of 100 /xm and 
a magnet with a transverse momentum kick of 200MeV/ c formed a spectrometer that mea- 
sured charged track momenta with resolution {ap/pY = (5x 10~^)^ + (1.4x 10~^p/(Ge"l//c))^. 
In-time energy deposits in the calorimeter were identified with a hardware cluster finder, and 
a fast energy sum for the entire calorimeter was used for triggering. In addition there were 
several annular photon veto counters surrounding the vacuum decay volume which enabled 
the rejection of events with decay products outside the acceptance region of the calorime- 
ter. Just behind the calorimeter was a 10 cm thick lead wall followed by a hodoscope bank 
of 45 overlapping counters designed to veto hadronic showers initiating in the calorimeter. 
The hadron veto was determined by an analog sum of all the individual signals from the 
counters. Following the hadron veto counter was a muon filter comprised of 3 m (20 interac- 
tion lengths) of steel. Muons traversing the filter steel were identified by a hodoscope bank 
consisting of 16 vertical, non-overlapping counters. 

In order to trigger on events with two muons we required the following elements: at least 
two hits in each of the trigger hodoscope banks; no hits in the photon veto counters; at 
least two non-adjacent hits in the muon-hodoscope bank; a sufficient number of hits in the 
drift chambers to be consistent with two tracks; no activity in the hadron veto counter; and 
a minimum energy sum of 6 GeV in the calorimeter. A two-track "minimum-bias" trigger 

4 



simultaneously collected Kl 7r~'~7r~7r° events which were used for normalization. This 
trigger was identical to the dimuon trigger but had no muon hodoscope, hadron hodoscope 
or minimum energy deposit requirements. The 7r"'"7r~7r° trigger was prescaled by a 

factor of 3600 relative to the n'^ trigger. 

Data were collected in the 1991 fixed target run during which 60 million dimuon triggers 
were recorded. Data reduction began by selecting events which had two good reconstructed 
tracks that projected to clusters in the calorimeter and one "photon cluster" not associated 
with tracks. These events were almost entirely Kl T^^fJ-^v^Jacc) or Kl 7r"'"7r~7r° 
decays. In the first case the pion is misidentified as a muon or the pion decayed, and "face is 
an "accidental" photon, i.e. not associated with the decay. For 7r"'"7r~7r° decays, 

both pions were misidentified as muons, and one of the photons from the decay of the 7r° 
escaped or overlapped a track cluster. The photon cluster was required to match a cluster 
found by the hardware cluster finder which had 30ns charge integration gate as opposed to 
the 100ns used to measure the energy deposited in each block. This suppressed events with 
out-of-time energy deposits. To further suppress contamination from accidental deposits, 
the photon cluster centroids were required to be at least 0.5 blocks from the beam holes 
in the calorimeter. Events with photon cluster energies below 8 GeV were rejected since 
accidental photons in background events typically have lower energy. This reduced the signal 
by 9% and rejected 42% of Kl T^^fJ-^v^jacc) events. Clusters not associated with tracks 
were required to have transverse profiles consistent with the shape expected from a photon. 
Track momenta were required to be above 7 GeV/c, greater than the 5 GeV/c threshold 
for muons to pass through 3 m of steel. Tracks with associated energy deposits above 3 
GeV in the calorimeter were rejected, thus suppressing events with hadron showers in the 
calorimeter. Monte Carlo studies of Kl n'^ events show that the photon clusters in 
background events are frequently near the beamholes in the calorimeter, and are therefore 
on average closer to tracks. We use this correlation to reject background events by requiring 
the separation of the photon cluster and the tracks at the calorimeter to be greater than 20 
cm. This cut reduces the signal by 11% and the Kl T^^fJ-^v^Jacc) background by 33%. 
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For events passing these cuts, the transverse momentum squared {Pf) versus the recon- 
structed invariant n'^ mass is shown in Fig. 1. There is a clear signal in the region near 
the Kl mass and low P^ . The invariant mass distribution for events with P'^ < 200MeV^/c^ 
is shown in Fig. 2. The region P^ < 200MeV^ / and n'^ invariant mass in the range 
0.482 < M^n^ < 0.514 GeV/c^ was chosen so as to optimize the statistical sensitivity to the 
branching ratio. For the Monte Carlo simulated data that passed the previous cuts, this 
signal box contains 77% of the n'^ events, and less than 0.5% of the background 

events. 

The large band at low /x"'"/x~7 mass in Fig. 1 is consistent with Kl 7r"'"7r~7r° decays, and 
it falls sharply with increasing mass, producing a negligible contribution under the signal 
peak. We also considered the radiative process Kl ir^fi^vj as a potential source of back- 
ground, however, it also falls rapidly with increasing mass and has negligible contribution 
at the kaon mass. The Kl T^^fJ-^v^jacc) background events produce a broad distribution 
in n'^ invariant mass and all of the background under the signal peak is consistent with 
having come from this process. All of the background processes considered have uniform 
P^ distributions, and we determined their combined mass-shape by using a region of data 
at 400 < P'} < 1000 MeV^c^ (Fig.2). The total background is estimated to be 10.5 ± 4.0 
events out of 207 candidate signal events. 

A similar analysis was performed on the minimum bias data to reconstruct Kl 7r"'"7r~7r° 
decays for normalization. The major difference is that now two photon clusters were required, 
and there was no 3 GeV track-cluster energy cut. The cut on the photon energy was relaxed 
to 7 GeV, thus increasing the yield. In addition, a cut of 3a (21 MeV) around the 7r° mass 
was made on the reconstructed invariant mass of the two photon clusters. 

High statistics Monte Carlo samples were generated for the normalization and signal 
modes. To properly account for activity in the detector not associated with the parent kaon 
decay, special "accidental" events were collected in parallel with the dimuon and minimum 
bias triggers. These accidental events sample the instantaneous rate of activity in the detec- 
tor throughout the run, and are overlaid on the simulated decays. This simulated data were 
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then passed through an identical analysis to that described above. Acceptances calculated 
by Monte Carlo for kaons decaying in the momentum range of 20-220 GeV/c and in a region 
90-160 m from the target were 1.95% for Kl M+M"7 and 1.95% for Kl 7r+7r-7r°. With 
196.5 signal events and a final sample of 20,919 7r"'"7r~7r° events, the branching ratio 

is B{Kl ^i+^i-j)/B{Kl 7r+7r-7r°) = (2.61 ± 0.19{stat)) x 10"^ 

The largest sources of systematic error in this ratio are due to differences between the 
dimuon and minimum bias triggers, the dominant contribution being from the hadron shower 
veto. Special muon calibration runs as well as minimum bias data were used to determine 
the absolute dimuon efficiency of the hadron veto to be (55.7 it 2.5)% [9]. The acceptance 
of the dimuon trigger hodoscope bank was (65.3 it 1.9)%, which includes contributions from 
geometric acceptance and trigger logic. The uncertainty in the determination of background 
was estimated to produce at most a 2.0% change in the branching fraction. After calibrating 
the lead glass there were some residual non-linearities, and this caused the branching fraction 
to have a small dependence on the photon cluster energy cut. The systematic uncertainty 
resulting from this cut was estimated to be 0.7%. Taking the branching fraction for 
7r"'"7r~7r° as 0.1238 it 0.0021 [11] and adding all errors in quadrature the branching fraction 
is B{Kl M+M"7) = (3.23 ± 0.23{stat) ± 0.19{sys)) x 10"^ This clearly estabhshes the 
decay, while agreeing with the only previous observation of one candidate event [12]. 

Our result corresponds to F^^^/F^^ = (5.66 it 0.59) X 10~^. Using Eq. 1 with f{x) = 1 
gives F^^^/F^^ = B{Kl /j.-'j)/ B{Kl 77) = 4.09 x 10"^ which is 2.7o- below our 

measured value. This is the first indication for the presence of a form factor in the decay 
Kl fi'^ The prediction of Ko [4], F^^^/F^^ = 7.45lo!i5 x 10~^ differs from our result 
by 2.90-. 

Assuming the model of Bergstrom, Masso, and Singer for the form factor, Eq. 1 can 
be integrated as a function of ax*- This is shown in Fig. 3, using identical expressions for 
/i(x) and fxix) as those of Ref. [6]. For this model our branching fraction corresponds to 
aK* = -0.018™. 

The spectrum of dimuon invariant masses is also sensitive to the form factor. In Fig. 4 



we show the dimuon mass distribution together with a Monte Carlo prediction for f{x) = 1. 
The background shape, normahzed to 10.5 events, was added to the Monte Carlo data. 
Again the background shape was determined using high-P^^ events. The assumption that 
the decay is dominated by the K"f"f* vertex leads directly to the dimuon invariant mass 
spectrum in Eq. 1, and the data in Fig. 4 follows this general shape, apart from the smaller 
contribution from f{x). Thus the data support the assumption of K"f"f* dominance. Also 
shown in Fig. 4 is the ratio of data to Monte Carlo, i.e. /{x)^. Using a maximum liklihood 
fit we obtain ax* = — O.lStoisB using the dimuon mass spectrum alone. Although less 
accurate, this is consistent with the value determined from the branching fraction. Together 
the combined best estimate is ax* = -0.028to:ni- This differs by 1.8a from the world 
average of ax* = —0.28 it 0.08 [11], using Kl e~^e~j. 

In conclusion, we have measured the branching fraction of the decay n'^ 
and the result is significantly above the prediction assuming a constant form factor. This 
supports the BMS model and is 2.9o" away from the model of Ko. The extracted form factor 
shows consistency with that from Kl e"'"e~7, at the 2a level. More data, particularly 
copious Kl e"'"e~7 with Kl n'^ will better test the viability of the BMS model. 
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computing and research divisions. One of us (Y.W.W) would like to acknowledge support 
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FIG. 1. 



vs — fj,'^ fj, 7 invariant mass distribution for events after cuts. 



11 



CM 



O 



60 



> 



50 



40 




LlJ 0.4 0.45 0.5 0.55 0.6 



FIG. 2. Distribution of 7 invariant mass after cuts. Histogram is data for events with 

< 200MeV^ /c'^ and points are data for events with 400 < < lOOOMeF^/c^, and normalized 




to the low P? data. 



12 



t 



CD 



I 

+ 

t 



x10 
120 



- 5 








OHL et. al. 
BARR et. ol. 



-1 



t 
BMS 
MODEL 



-0.5 0.5 1 

FIG. 3. This result for the branching ratio of /x"'"/x~7, normalized to 77 (large 

black dot) and the extracted range of the parameter a^* from the model of Bergstrom, Masso, 
Singer [1]. Also shown (horizontal error bars) are the extracted values of a^* from Refs. [6] and 
[5] for the case of e~^e~j. 
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FIG. 4. Top: Dimuon mass distribution, points are data, histogram is a Monte Carlo simulation 
with /(cc) = 1, dashed line is the estimated background, dotted line is the acceptance. Bottom: 
the ratio of data to Monte Carlo from the top plot (i.e. our measured values of /(cc)^). The solid 
line is the result from a maximum liklihood fit. For illustration we also show predictions assuming 
ax* = — 1 (dashed) and ax* = +1 (dotted). For these plots the normalization was fixed by the 
number of events. 



14 



